flD-fll59  221  COLOR  POST-PROCESSOR  FOR  THE  EPIC-2  COHPUTER  PROGRAM  1/1 

ARHV  BALLISTIC  RESEARCH  LAB  ABERDEEN  PROVING  GROUND 
HD  T  M  SHERRICJf  JUN  85  BRL-5P-48 


UNCLASSIFIED 


F/G  9/2 


NL 


.  I-'W- 


'--T  rj  <rrjrjLi  ‘j  y j.yj  yj’ir  ]  •.- 


MICROCOPY  RESOLUTION  TEST  CHART 

national  bureau  OE  STAN0AR0S-1963-A 


•yy  vj  iTjy  '  yj  y  ji:  \m<’m  ■■■jJvmj'uiiMJM 


AD-A159  321 


AD 


US  ARMY 
MATERIEL 
COMMAND 


SPECIAL  PUBLICATION  BRL-SP-48 


COLOR  POST-PROCESSOR  FOR  THE  EPIC-2 

COMPUTER  PROGRAM 

Thomas  M.  Sherrick 

June  1985 

APPROVED  FOR  PUBUC  RELEASE;  DISTRIBUTION  UNUMITED. 

US  ARMY  BALLISTIC  RESEARCH  LABORATORY 

ABERDEEN  PROVING  GROUND,  MARYLAND 


DTIC 

ELECTEj 

SEP  2  3  1985 


"Original  containn  color 
plates:  All  7 riC  reproduct¬ 
ions  will  be  in  black  and 
White" 


03  0? 


85  9 


Destroy  this  report  when  it  is  no  longer  needed. 
Do  not  return  it  to  the  originator. 


Additional  copies  of  this  report  may  be  obtained 
from  the  National  Technical  Information  Service, 

U.  S.  Department  of  Commerce,  Springfield,  Virginia 


22161. 


The  findings  in  this  report  are  not  to  be  construed  as  an  official 
Department  of  the  Army  position,  unless  so  designated  by  other 
authorized  documents. 


The  use  of  trade  names  or  manufacturers'  names  in  this  report 
does  not  constitute  indorsement  of  any  commercial  product. 


IN 


security  classification  of  this  page  (When  Dmtm  Entered) 


REPORT  DOCUMENTATION  PAGE 


t.  REPORT  NUMBER 


12.  GOVT  ACCESSION  NO. 


Special  Publication  BRL-SP-48 


4.  TITLE  cand  Subtitle; 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3.  RECIPIENT'S  CATALOG  NUMBER 


5.  TYPE  OF  REPORT  &  PERIOD  COVERED 


Color  Post-Processor  for  the  EPIC-2  Computer 
Program 


7.  AUTHORfe; 


Final 


6.  PERFORMING  ORG.  REPORT  NUMBER 


8.  CONTRACT  OR  GRANT  NUMBERCa) 


Thomas  M.  Sherrick 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

US  Army  Ballistic  Research  Laboratory 
ATTN:  AMXBR-TBD 

Aberdeen  Proving  Ground,  MD  21005-5066 


II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

US  Army  Ballistic  Research  Laboratory 
ATTN:  AMXBR-OD-ST 

Aberdeen  Provine  Ground.  MD  21005-5066  _ 


4.  monitoring  agency  name  a  A00RESS('I/  dliferent  from  Controiting  Office) 


10.  PROGRAM  ELEMENT,  PROJECT,  TASK 
AREA  &  WORK  UNIT  NUMBERS 


IL162618AH80 


12.  report  date 
June  1985 

13.  NUMBER  OF  PAGES 

.50 


IS.  SECURITY  CLASS,  (ol  thia  report; 


IS.  DISTRIBUTION  STATEMENT  Co/  tbia  Report; 


UNCLASSIFIED 


ISa.  DECLASSIFICATION/ downgrading 
SCHEDULE 


Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  statement  (ol  tba  abatraet  antarad  In  Block  20.  II  dlltormt  from  Report) 


18.  SUPPLEMENTARY  NOTES 


"Original  contains  color 
plates:  All  DTIC  reproduct¬ 
ions  will  be  in  black  and 
white" 


19.  KEY  WORDS  (Continue  on  reveree  eide  if  neceBemry  end  identify  by  block  number) 


state  Plots 
Color  Graphics 
Post-Processor 
Contours 


Shading 

Velocity  Vectors 


2(L  ABSTRACT  (*Ctet£BU«  em  rer^ee  sM*  ff  rmeeeemty  mod.  idenlity  by  block  number) 

The  Color  Post-Processor  for  the  Epic-2  Computer  Program  is  an  extensive 
modification  of  the  original  Post-Processor  for  the  Epic-2  Computer  Program. 
This  report  documents  modifications  made  to  the  State  Plots  Program  in  the 
Post-Processor  and  the  implementation  of  color  graphics  into  the  Post-Processoi 
The  State  Plots  Program  represents  the  voluminous  amount  of  numerical  data 
generated  by  the  Epic-2  Computer  Program  as  Deformed  Geometry,  Stress  Field, 
Pressure  Field,  Strain  Field  and  Velocity  Vector  Plots.  The  color  graphics 
incorporated  into  the  Post-Processor  illustrate  shading  an 


DD  t  JAN  7>  1473  eOfTIOM  or  f  NOV  6S  IS  OBSOLETE  UNCLASSIFIED 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  fl*»>««  Dmim  Entorod) 


UNCLASSIFIED _ 

security  classification  of  this  PAOEfini«n  Bnl«f«f> _ 

20.  Continued... 

the  respective  plots  in  the  State  Plots  Program.  Detailed  instructions  for 
using  the  Color  Post-Processor  are  included. 


UNCLASSIFIED 


security  classification  of  this  PAGEflFhan  Data  Enfarad; 


TABLE  OF  CONTENTS 


LIST  OF  FIGURES  .... 

I .  INTRODUCTION  . 

II.  THE  POST-PROCESSOR  .  .  . 

III.  CODE  MODIFICATIONS  .  . 

IV.  PROGRAM  USER  INSTRUCTIONS 

V.  SUMMARY  OF  INPUT  DATA  . 

VI .  SUMMARY  . 

DISTRIBUTION  LIST  .  .  . 


Accession  For 


NTIS  GRA4I  1^ 

DTIC  TAB  □ 

Unannovinced 
Justification _ 


i  By - 

!  Distribution/ 


I  Availability  CoJos 
jAvail  and/or 
I  lat  I  Special 


LIST  OF  FIGURES 


rWJ  IT  ■  B“J  W  :•  1 


Figure 


Pressure  Contours  for  Material  1  only  (Projectile)  with 

IDEM  =  0 . 10 

Pressure  Contours  for  Material  2  only  (Target)  with  IDEM  =  5  ...  10 


Velocity  Vectors  with  IDEM  =  0 


Velocity  Vectors  with  IDEM  =  5 


Velocity  Vector  Plot . 11 


Axial  Load  Plot 


Stress  Color  Shaded  Plot 


Pressure  Color  Shaded  Plot  . 

Strain  Color  Shaded  Plot  . 

Stress  Color  Contour  Plot  . 


Pressure  Color  Contour  Plot 


Strain  Color  Contour  Plot 


I .  INTRODUCTION 


This  report  documents  enhancements  made  to  the  EPIC-2  post-processor  and 
the  implementation  of  color  into  the  post-processor.  Reference  1  describes 
the  EPIC-2  program  and  contains  a  user's  manual  for  the  program  and  post¬ 
processor.  This  report  can  be  utilized  as  documentation  for  the  original 
EPIC-2  post-processor,  STPLOTPL,  by  merely  neglecting  the  information  per¬ 
taining  to  color  graphics. 

Plotting  information  for  the  post-processor  is  written  on  two  tapes  by 
the  EPIC- 2  program  at  user-specified  intervals.  One  tape  contains  all  the  flow 
field  variables  at  preselected  times.  The  other  contains  time  histories  of 
flow  field  variable  at  preselected  times  and  locations. 

As  described  in  [1],  EPIC-2  is  a  two-dimensional  computer  program  for 
impact  and  explosive  detonation  problems.  The  capability  to  treat  spinning 
projectiles  is  also  included.  The  numerical  technique  is  based  on  a  Lagrangian 
finite  element  formulation  in  which  the  equations  of  motion  are  integrated 
directly  rather  than  through  the  traditional  stiffness  matrix  approach.  Tri¬ 
angular  elements  are  formulated  for  large  strains  and  displacements,  and  non¬ 
linear  material  strength  compressibility  effects  are  included  to  account  for 
elastic  plastic  flow  and  wave  propagation.  ,  ,  ^ 

II.  THE  POST- PROCESSOR 

The  EPIC-2  post-processor  is  a  plotting  program  which  reads  data  from 
previously  generated  tapes  in  an  EPIC-2  calculation.  The  post-processor  can 
produce  state  plots  and  time  plots. 

The  state  plots  are  broken  down  into  five  types: 

1.  Deformed  geometry  plots  -  plots  composed  of  triangular  or  triangle 
in  quadilateral  elements. 

2.  Stress  field  plots  -  plots  the  lines  of  constant  equivalent  stress 
given  by  the  following  equation; 

5  =  y  3/2  ♦  s|  .  .  3  .  T|g)  .  0) 

3.  Pressure  field  plots'-  plots  net  pressure,  that  is,  hydrostatic 
pressure  plus  artificial  viscosity. 

4.  Strain  field  plots  -  plots  equivalent  plastic  strain  given  by  the 
following  equation; 


-  t  +  At 
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1 

Johnson^  G.  E.,  "EPIC-2j  A  Computer  Program  for  Elastia-Plastta  Impact 
Computations  in  Tuo  Dimensions  Plus  Spin,"  ARBRL-CR-00373y  June  1978 
(ADA058786). 


where 


The  post-processing  program,  as  received,  proved  to  be  very  inefficient. 
Changes  were  made  to  various  sections  of  the  state  plot  portion  of  the  post¬ 
processor.  At  the  same  time,  color  capability  was  also  added  as  the  result  of 
acquisition  of  Tektronix  4027  color  graphics  terminals.  These  changes  are 
described  below.  Modification  of  the  time  plot  portion  of  the  post-processor 
will  be  described  in  a  future  report. 


Ill,  CODE  MODIFICATIONS 


Subroutine  EDGE  proved  to  be  the  most  time-consuming  part  of  the  original 
EPIC-2  post-processor.  This  subroutine  plots  the  outline  of  the  figure  being 
plotted.  It  was  modified  by  incorporating  a  sort/merge  routine  into  the  edge 
subroutine,  thus  reducing  the  computing  costs  by  33%  to  50%.  This  subroutine 
is  used  in  a  large  portion  of  the  post-processor,  that  is,  it  is  used  for 
deformed  geometry  plots,  pressure  field  plots,  stress  field  plots  and  velocity 
field  plots. 

A  unit  conversion  subroutine  was  added  to  the  post-processor.  In  this 
subroutine,  SIUNIT,  variables  are  converted  from  Engineering  Units  to  Standard 
International  Units.  There  is  also  an  entry  point  in  the  subroutine,  EGUNIT, 
which  converts  Standard  International  Units  to  Engineering  Units.  This  sub¬ 
routine  is  invoked  by  setting  ICONV  to  some  integer  value  on  the  second  input 
card.  This  card  is  explained  in  more  detail  in  the  input  instructions  section 
of  this  report. 

Material  selection  and  density  parameters  were  added  to  the  pressure,  stress, 
strain,  shear  and  velocity  field  plots.  The  material  selection  parameter 
allows  for  the  plotting  of  contours  or  velocity  vectors  for  a  user  specified 
material  when  working  with  a  multi-material  calculation.  For  example,  for  the 
case  of  a  projectile  penetrating  a  target  of  different  material,  material 
selection  permits  the  contours  to  be  displayed  in  either  the  projectile  or  the 
target  or  both.  The  density  parameter,  IDEN,  determines  the  number  of  symbols 
to  be  plotted  on  the  contour  itself,  or  the  number  of  velocity  vectors  to  be 
plotted.  The  IDEN  parameter  is  explained  more  fully  in  the  input  instructions 
section.  These  features  are  illustrated  in  Figures  1-4. 

A  feature  incorporated  into  the  velocity  vector  subroutine  is  the  addition 
of  arrowheads  to  the  velocity  vectors  in  order  to  show  direction.  Figure  5 
is  a  typical  example  of  a  velocity  vector  plot. 

In  each  of  the  subroutines  PPLOT,  SPLOT,  EPLOT,  and  CSHEAR,  there  is  a 
process  which  performs  four  tasks; 

1.  Checks  the  elements  to  see  if  they  are  beyond  the  extreme  coordi¬ 
nates  or  completely  failed. 

2.  Determines  if  the  pressure,  stress,  shear  or  strain  contour  is 
present  between  each  pair  of  nodes  in  the  element. 

3.  Calculates  the  position  of  the  pressure,  stress,  shear  or  strain 
contour  intersection  on  a  line  between  a  pair  of  nodes. 

A.  Fills  arrays  RP  and  ZP  for  drawing  the  pressure,  stress,  shear 
or  strain  contour. 

Basically,  this  algorithm  loops  through  the  elements  and  determines  the  nodes 
associated  with  each  element.  The  algorithm  then  averages  the  pressures  of 
the  elements  surrounding  each  node  and  linearly  interpolates  between  the  nodes 
to  find  where  the  contour  intersects  the  element  boundary.  The  process  has 
now  been  incorporated  into  a  subroutine  called  CONPLOT,  thus  making  the 
program  more  maintainable  and  less  time  consuming. 
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FIGURE  1--Pressure  Contours  for  ma¬ 
terial  1  only  (Projectile)  with  IDEN  =  0 


FIGURE  3--Velocity  Vectors  with  IDEN  =  0 
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FIGURE  2--Pressure  Contours  for  Ma¬ 
terial  2  only  (Target)  with  IDEN  =  5 


FIGURE  4--Velocity  Vectors  with  IDEN  =  5 


FIGURE  5--Velocity  Vector  Plot 


FIGURE  6--Axial  Load  Plot 
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FIGURE  7--Stress  Color  Shaded  Plot 


FIGURE  8--Pressure  Color  Shaded  Plot 


In  addition,  this  subroutine  calculates  the  maximum  value  of  stress,  strain, 
pressure  and  shear  using  the  standard  FORTRAN  intrinsic,  min-max  function 
routines.  After  finding  the  maximum  value,  five  contour  lines  are  set  up 
plotting  98%,  80%,  60%,  40%,  and  20%,  respectively,  of  the  maximum  constant 
stress,  pressure,  strain  or  shear. 

Subroutine  LOAD  was  also  added  to  the  post-processor.  This  subroutine 
calculates  and  plots  the  axial  load  in  the  projectile.  This  is  accomplished 
by  calculating  the  force  in  the  rings  of  elements  of  the  projectile,  multiply¬ 
ing  pressure  times  area.  This  can  be  obtained  only  along  the  centerline  and 
is  only  meaningful  for  normal  impacts.  Axial  load  can  be  an  especially  useful 
feature  for  situations  involving  long  rod  impacts.  An  example  of  axial  load 
plots  for  an  L/D  =  25  steel  projectile  striking  a  rigid  barrier  is  shown  in 
Figure  6. 

Principal  stress  and  maximum  shear  contours  may  also  be  plotted  by  the 
EPIC-2  post-processor.  This  is  a  color  contour  representation.  The  following 
equations  are  evaluated  in  the  READ  2  subroutine  in  the  post-processor; 

1.  Maximum  shear; 


2.  Principal  stresses; 


<7 

r  ,z 


C5) 


STRESSES  ARE  FOR  THE  (r-z)  PLANE. 

To  implement  the  shear  capability,  another  block  of  storage,  called  SHEAR, 
was  added  to  the  common  block  /ELEM/.  Care  must  be  taken  so  as  not  to  exceed 
the  common  block  length.  This  is  dependent  upon  the  memory  capacity  of  the 
machine  running  the  program. 


In  order  to  incorporate  color  into  the  velocity  vector  plotting  capability 
of  the  EPIC-2  post-processor,  it  was  necessary  to  make  some  changes  to  the 
subroutine  VPLOT.  The  velocity  vectors  are  scaled  so  as  to  give  a  vector  one 
inch  long  to  represent  the  maximum  velocity  in  the  grid  for  both  raster  and 
vector  terminals.  On  the  Tektronix  4027  terminal  one  inch  is  equivalent  to  64 
rasters.  The  scaling  is  done  in  the  following  manner; 


XT  =  R(I} 
vector. 


XT,  YT  represent  the  coordinates  of  the  tail  of  the  velocity 
R(I)  IS  THE  RADIAL  COORDINATE  OF  THE  NODE. 


YT  =  Zfl)  Z(I)  IS  THE  AXIAL  COORDINATE  OF  THE  NODE. 


Call  GTSCAL  (SFX,  SFY)  -  this  routine  obtains  the  proper  scale  factor  and 
converts  the  R,  Z  values  to  user  units  per  raster. 
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VI .  SUMMARY 


The  EPIC-2  post-processor  color  program  is  currently  available  at  the 
Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground,  Maryland.  It  is  set  up 
to  run  with  the  EPIC- 2  computer  program  on  the  Cyber  7600.  The  output  is  dis¬ 
played  on  the  Tektronix  4027  color  graphics  terminal  and  hard  copy  is  obtained 
from  a  Matrix  color  hard  copy  recording  system,  model  4007. 

Future  plans  include  modification  of  the  time  plots  program  and  implementing 
color  into  this  capability. 


V.  SUMMARY  OF  INPUT  DATA 

BLANK  CARD 

State  plot  identification  card.  One  required  for  each  type  of  plot  desired 
(215,  211,  3X,  4F10.0). 

TYPE  CYCLE  ICAL  ICONV  ZMAX  ZMIN  RMAX  RMIN 
Deformed  geometry  plot  card.  Type  =  1  C3A6,  2X,  5I5X 
TITLE  EDGE  IE  ID  IF  IM 

Stress,  pressure,  strain  or  shear  plot  card.  Type  =  2,  3,  4,  or  6.  Two  cards 
required.  (Contour  plots) 

Card  #1  (3A6,  2X,  215,  5F10.0) 

TITLE  EDGE  NLINE  PARAM(l)  PARAM(2)  PARAM(3)  PARAM(4)  PARAM(5) 

Card  #2  (215) 

I DEN  IMAT 

Stress,  pressure,  strain  or  shear  plot  card.  Type  =  7,  8,  9  or  10.  Two  cards 
are  required.  (Shaded  plots) 

Card  #1  (3A6) 

TITLE 

Card  #2  (5X,  15) 

IMAT 

Velocity  vector  plot  card.  Type  =  5  (3A6,  2X,  215,  FIO.O,  15) 

TITLE  EDGE  NSTRT  VSCALE  I DEN 

Axial  load  plot  card.  Type  =  11.  Only  state  plot  identification  card  required. 
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Card  #1  (3A6,  2X,  215,  5F10.0) 

TITLE  (3A6)  -  Title  to  be  written  on  the  plot.  Presently  there  is  no  need 
for  the  remainder  of  the  fields  on  this  card.  The  EDGE  sub¬ 
routine  is  not  required  because  the  shaded  plots  generate  their 
own  edge.  Also,  the  NLINE  value  and  PARAM(I)  values  are  unnec¬ 
essary,  because  the  shading  subroutines  generate  their  own  de¬ 
grees  of  shading.  The  format  remains  the  same  for  further  de¬ 
velopment. 

Card  #2  (215) 


The  IDEN  parameter,  the  first  15  field  on  this  card,  is  not  needed  since  no 
symbols  are  plotted  on  the  shaded  plots. 

IMAT  (15)  -  Plot  shaded  areas  for  IMAT,  if  IMAT  =  0,  all  materials  are 

plotted. 

Present  only  if  TYPE  =  5 


Velocity  vector  plot  card  (3A6,  2X,  215,  FIO.O,  15) 


This  type  of  plot  generates  velocity  vectors  based  on  the  position,  velocity  and 
direction  of  the  velocity  of  the  nodes. 

TITLE  (3A6)  -  The  title  which  is  written  on  the  plot. 


EDGE  (15)  - 


Number  to  indicate  type  of  outline  to  be  plotted,  same  as  those 
described  for  the  geometry  plot. 


NSTRT  (15)  - 


The  node  number  to  start  plotting  the  velocity  vectors,  if 
NSTRT  =  0  all  nodal  velocity  vectors  are  plotted. 


VSCALE  (FIO.O)  -  That  velocity  vector  which  produces  a  velocity  vector  of  length 
1.0  inch. 

IDEN  (15)  -  Density  control  parameter  for  velocity  vectors,  for  IDEN  =  0 

all  velocity  vectors  are  plotted.  This  is  based  on  the  same 
type  modulo  function  as  the  other  IDEN  parameters. 
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EDGE  =  -1  plots  an  outline  of  the  mirror  image  of  the  actual  geometry  for 
axisymmetric  cases. 

EDGE  =  2  plots  an  outline  of  both  actual  geometry  and  the  mirror  of  actual 

geometry. 

EDGE  =  3  plots  an  outline  of  the  entire  image  without  showing  the  elements. 

IE (15)  -  =1  prints  E  in  those  elements  which  are  in  the  elastic  range. 

ID(I5)  -  =1  prints  P  in  those  elements  which  are  in  the  plastic  flow 

range. 

IF(I5)  -  =1  prints  F  in  those  elements  which  have  failed. 

IM(15)  -  =1  prints  material  number  inside  elements.  Effective  only  if 

IE  =  ID  =  IF  =  0  or  are  left  blank. 

Present  only  if  TYPE  =  2,  3,  4,  and  6.  These  types  give  plots  of  constant 
equivalent  stress,  net  pressure,  equivalent  plastic  strain  or  shear  contours. 

There  are  two  cards  required  for  each  of  these  plots. 

Card  #1  (3A6,  2X,  215,  5F10.0) 

TITLE  (3A6)  -  Title  to  be  written  on  the  plot. 

EDGE  (15)  -  Number  to  indicate  type  of  outline  to  be  plotted.  (Same  as 

those  described  for  the  geometry  plot.) 

NLINE  (15)  -  The  number  of  contour  lines  to  be  plotted,  ^fust  be  1,  2,  3,  4 
or  5. 

PARAM  (I)  -  Specific  contour  values  of  stress,  net  pressure,  equivalent 

1=1,  NLINE  plastic  strains  or  shear.  There  is  an  algorithm  which  will 

(FIO.O)  automatically  calculate  98%,  80%,  60%,  40%,  and  20%,  respec¬ 

tively,  of  the  maximum  value  of  the  chosen  variable.  If  this 
capability  is  chosen,  then  0.0  is  entered  in  the  first  (FIO.O) 
field  of  the  five  and  the  remaining  four  fields  are  left  blank. 

Card  #2  (215) 

IDEN  (15)  -  Density  of  contour  label  symbols.  If  (MOD  (Node  Number/ Element 

Number,  IDEN))  is  0,  a  label  symbol  will  be  plotted. 

IMAT  (15)  -  Plot  contours  only  for  IMAT,  material  selection  parameter,  if 

IMAT  =  0,  all  materials  are  plotted. 

Present  only  if  TYPE  =  7,  8,  9,  and  10.  These  plots  give  shaded  areas  of 
Stress,  Pressure,  Strain,  or  Shear.  The  areas  are  determined  by  the  respective 
subroutines  as  mentioned  earlier  in  the  report. 

Again,  two  cards  are  required  for  each  of  these  plots. 


TYPE  =  7  plots  a  shaded  stress  field 

TYPE  =  8  plots  a  shaded  pressure  field 

TYPE  =  9  plots  a  shaded  strain  field 

TYPE  =  10  plots  a  shaded  shear  field 

TYPE  =  11  plots  the  axial  load  distribution  in  the  projectile. 

CYCLE  (15)  -  The  cycle  number  which  is  to  be  plotted.  The  printed  output  of 
the  main  program  gives  the  cycle  numbers  of  the  data  written  on 
tape  2. 

ICAL  (II)  -  This  is  a  CALCOMP  option  and,  in  this  instance,  will  always  be 

0. 

ICONV  (II)  -  If  ICONV  =  0,  the  output  (coordinates,  velocities,  stresses  and 
pressures)  will  be  in  SI  units.  If  ICONV  is  greater  than  0,  the 
output  will  be  in  Engineering  units. 

ZMAX  (FIO.O)  -  The  maximum  axial  Z  coordinate  on  the  plot.  The  Z  axis  is  10 
inches  long  and  is  divided  into  10  equal  increments  from  ZMAX 
to  ZMIN.  ZMAX  represents  the  maximum  Z  coordinate  of  the  phys¬ 
ical  problem  analyzed  by  EPIC- 2.  Ignored  for  TYPE  =  11  plots. 

ZMIN  (FIO.O)  -  The  minimum  Z  coordinate  on  the  plot.  Ignored  for  TYPE  =  11 
plots. 

RMAX  (FIO.O)  -  The  maximum  radial  R-coordinate  on  the  plot.  The  R  axis  (RMAX 

to  RMIN)  can  be  any  length  and  the  scale  is  equal  to  that  of  the 
Z  axis.  RMAX  represents  the  maximum  R-coordinate  of  the  physi¬ 
cal  problem  analyzed  with  EPIC-2. 

RMIN  (FIO.O)  -  The  minimum  R  coordinate  on  the  plot.  Ignored  for  TYPE  =  11 
plots. 

The  next  card  required  describes,  with  more  detail,  the  type  of  plot  chosen 

by  the  first  card.  The  format  is  dependent  upon  the  type  of  plot  chosen. 

Present  only  if  TYPE  =  1 

Deformed  geometry  plot  card  (3A6,  2X,  515) 

TITLE  (3A6)  -  The  title  which  is  written  on  the  plot.  The  time  and  cycle  num¬ 
ber  are  also  written  on  the  plot  using  data  previously  read  by 
the  post-processor  subroutines. 

EDGE  (15)  -  A  number  to  specify  if  an  outline  of  the  geometry  should  be 

plotted. 

EDGE  =  0  plots  no  outline. 

EDGE  =  1  plots  an  outline  of  the  actual  geometry. 


FIGURE  9--Strain  Color  Shaded  Plot 


FIGURE  C-2 — Pressure  Color  Contour  Plot 


FIGURE  C-1--Stress  Color  Contour  Plot 


FIGURE  C-3 — Strain  Color  Contour  Plot 
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set  up,  the  ranges  in  the  algorithm  are; 


IF 

(FAC  . 

GE. 

.5)  ICOLOR  =  1 

IF 

(.36  . 

LE. 

FAC  .AND.  FAC  .LT. 

.5)  ICOLOR  =  2 

IF 

(.25  . 

LE. 

FAC  .AND.  FAC  .LT. 

.36) 

ICOLOR  =  3 

IF 

(.14  . 

LE. 

FAC  .AND.  FAC  .LT. 

.25) 

ICOLOR  =  4 

IF 

(.036 

.LE, 

,  FAC  .AND.  FAC  .LT. 

.14) 

ICOLOR  = 

IF 

(.002 

.LE, 

.  FAC  .AND.  FAC  .LT. 

.036) 

ICOLOR  = 

IF 

(FAC  . 

,LT. 

.002)  ICOLOR  =  0 

where  FAC  is  the  normalized  value. 

2.  In  the  second  algorithm,  threshold  values  of  pressure,  stress  or  shear 
are  used  as  the  limiting  factor.  For  example,  consider  the  handling  of  com¬ 
pression  and  tension  with  respect  to  a  shaded  pressure  PLOT,  one  color  can  rep¬ 
resent  all  positive  pressures,  compression,  another  could  represent  all  negative 
pressures,  tension,  finally  a  third  could  represent  those  pressures  which  are 
essentially  zero.  This  algorithm  uses  the  actual  values  of  the  pressure  from 
the  tape  generated  by  the  EPIC- 2  program.  These  values  are  not  normalized. 

The  shaded  plot  capability  is  shown  in  figures  7-9. 

IV.  PROGRAM  USER  INSTRUCTIONS 

The  input  required  for  the  color  processor  is  similar  to  that  of  the  old 
post-processor.  However,  for  the  sake  of  completeness,  it  is  worth  reviewing 
all  the  input.  Some  portions  of  the  following  are  taken  directly  from  the  old 
post-processor  description  found  in  [1]. 

The  first  card  is  blank.  This  card  is  used  for  CALCOMP  information  and  is 
not  used  by  the  color  processor.  However,  it  is  required. 

The  next  card  is  a  state  plot  identification  card.  The  format  is  (215,  211, 
3X,  4F10.0)  representing: 

TYPE  (15)  -  a  number  indicating  the  t)'pe  of  plot  desired. 

TYPE  =  1  plots  a  geometry  plot  i.e.  element  connectivity 

TYPE  =  2  plots  stress  field  contours 

TYPE  =  3  plots  pressure  field  contours 

TYPE  =  4  plots  strain  field  contours 

TYPE  =  5  plots  a  velocity  field 

TYPE  =  6  plots  shear  field  contours 


15.  SUBROUTINE  DELAY  CIMSEC)  -  causes  the  4027  to  wait  at  least  /IMSEC/ 
milliseconds. 

16.  SUBROUTINE  BELL27  -  sounds  the  4027  bell. 

17.  SUBROUTINE  BAUD  (ITR,  IR)  -  sets  the  transmit  and  receive  Baud. 

18.  SUBROUTINE  RMAP  CICOLOR,  IHUE,  ILIGHT,  ISAT)  -  changes  color  /ICOLOR/ 
by  relative  values  of  hue,  lightness  and  saturation. 

19.  SUBROUTINE  RPOL  (IX,  lY,  NPT.,j  -  creates  a  polygon  with  vertices  given 
in  relative  coordinates. 


20.  SUBROUTINE  RVE  (IX,  lY,  NPTS)  -  draws  a  vector  in  relative  coordinates. 

21.  SUBROUTINE  JUMP  (IROW,  ICOLM)  -  moves  the  workspace  Alpha  cursor  to  row 
/IRON/  and  columns  /ICOLM/  of  the  workspace. 

22.  SUBROUTINE  GRA  (IROWl,  IROW2,  ICOLMl,  IC0LM2)  -  creates  a  graphic  area 
in  the  workspace  containing  rows  /IROWl/  through  /IR0W2/  of  columns  /ICOLMl/ 
through  /IC0LM2/. 

23.  SUBROUTINE  ANWORD  (NCHAR,  IWORD)  -  this  subroutine  outputs  up  to  2 
alphanumeric  words  to  the  CRT. 

24.  SUBROUTINE  DISPNO  (FNUMB,  IDEC)  -  displays  the  values  for  the  contours. 

25.  SUBROUTINE  SHRINK  -  shrinks  graph  address  so  that  the  4027  can  emulate 
a  4010  terminal. 


26.  SUBROUTINE  COPY27  (ITYPE)  -  this  subroutine  directs  the  4027  to  copy  a 
host  file  to  the  4924  tape  drive  with  simultaneous  display,  ITYPE  =  0,  or  with¬ 
out  display,  ITYPE  =  1. 

A  major  contribution  to  the  color  post-processor  is  the  creation  of  the  shaded 
state  plots  for  stress,  pressure,  strain  and  shear  fields.  The  subroutines  are 
basically  similar,  hence,  the  description  of  one  subroutine  is  sufficient  to 
demonstrate  the  process  involved. 


The  subroutine  begins  by  setting  the  eight  basic  colors,  which  can  be 
plotted  simultaneously  occupying  distinct  regions,  to  the  values  of  whichever 
colors  are  desired.  The  4027  has  eight  basic  colors,  CO  through  C7  which  are 
white,  red,  green,  blue,  yellow,  cyan,  magenta  and  black,  respectively.  By 
using  subroutine  MAP  these  colors  can  be  changed  to  any  of  the  128  colors  capa¬ 
ble  of  being  produced  by  the  4027  color  terminal.  The  maximum  value  of  stress, 
pressure,  strain  or  shear  is  then  calculated  using  the  standard  FORTRAN  intrin 
sic,  min-max  function  routines.  Finally,  one  of  two  algorithms  is  employed  to 
determine  the  anwunt  of  shading  desired.  The  two  algorithms  are: 


1.  The  first  algorithm  normalizes  all  the  values  for  pressure,  stress, 
strain  or  shear  by  dividing  by  the  maximum  value  obtained  from  the  maximizing 
function.  Ranges  of  values  are  then  set  up  which  can  produce  as  many  or  as 
few  colors  as  desired.  For  example,  the  way  the  program  is  currently 


XH  =  XT  +  (RDOT(I)/VSCALE)  *  64.  *  SFX  »  ROOT  Is  radial  velocity  of  the  node. 

YH  =  YT  +  (ZDOT(I)/VSCALE)  *  64.  *  SPY  ,  ZDOT  is  axial  velocity  of  the  node. 

VSCALE  is  the  velocity  scaling  factor. 

XH,  YH  represent  the  coordinates  of  the  head  of  the  velocity  vector; 
these  values  are  scaled  to  fit  the  4027  color  terminal. 

After  scaling  is  completed,  a  test  is  performed  to  check  if  the  vector  length 
is  less  than  .05  inches.  If  so,  it  is  not  plotted.  Finally  an  arrowhead  is 
placed  on  the  vector  denoting  its  direction  and  the  scale  for  the  vector  is 
displayed  on  the  screen. 

In  implementing  color  into  the  post-processor,  a  color  library  was  created. 
It  is  composed  of  subroutines  completed  with  PLOTIO  and  the  Tektronix  4027  color 
graphics  terminal.  The  library  subroutines  are  written  in  FORTRAN  and  allow 
the  user  control  of  the  4027  color  terminal  with  simple  CALL  commands.  The 
library  contains  the  following  subroutines: 

1.  SUBROUTINE  DRWLIN  -  draws  a  line. 

2.  SUBROUTINE  DRWARR  -  draws  lines. 

3.  SUBROUTINE  CIR(IR)  -  draws  a  circle  with  radius  IR. 

4.  SUBROUTINE  COL(ICOLOR)  -  sets  the  current  vector  and  polygon  to  ICOLOR. 

5.  SUBROUTINE  ERA  -  erases  the  workspace,  monitor  or  graphics  region. 

6.  SUBROUTINE  ERASEG  (ICOLOR)  -  erases  the  contents  of  the  graphic  area 
with  ICOLOR. 

7.  SUBROUTINE  LIN(ITYPE)  -  select  line  type,  i.e.,  solid,  dashed,  etc. 

8.  SUBROUTINE  MAP  (ICOLOR,  IHUE,  ILIGHT,  ISAT)  -  sets  color,  ICOLOR,  to 
hue,  lightness,  and  saturation  levels  contained  in  /IHUE/,  /ILIGHT/,  and  /ISAT/, 
respectively. 

9.  SUBROUTINE  MIX  (ICOLOR,  IRED,  IGREEN,  I  BLUE)  -  sets  color,  ICOLOR,  to  a 
given  color  value  expressed  as  percentages  of  red,  green,  and  blue. 

10.  SUBROUTINE  PIE  (IR,  ISTART,  ISTOP)  -  draws  a  pie  with  a  radius,  IR, 
fills  in  the  pie  area  between  /ISTART/  and  /ISTOP/,  which  are  expressed  in 
in  degrees. 

11.  SUBROUTINE  POL  (IX,  lY,  NPTS)  -  color  fills  a  polygon  with  /NPTS/  ver¬ 
tices. 

12.  SUBROUTINE  VEC  (IX,  lY,  NPTS)  -  draws  a  vector  in  the  graphics  region, 

13.  SUBROUTINE  ENA  -  places  the  4027  in  the  Gin  (GRAPHIC  INPUT)  mode, 

14.  SUBROUTINE  DISA  -  causes  the  4027  to  exit  the  Gin  mode. 
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Department  of  Physics 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the 
reports  it  publishes.  Your  comments/answers  to  the  items/questions  below  will 
aid  us  in  our  efforts. 

1 .  BRL  Report  Number _ Date  of  Report _ 

2.  Date  Report  Received _ 

3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or 

other  area  of  interest  for  which  the  report  will  be  used.) _ 


4.  How  specifically,  is  the  report  being  used?  (Information  source,  design 
data,  procedure,  source  of  ideas,  etc.) _ 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far 
as  man-hours  or  dollars  saved,  operating  costs  avoided  or  efficiencies  achieved, 
etc?  If  so,  please  elaborate. _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future 
reports?  (Indicate  changes  to  organization,  technical  content,  format,  etc.) 
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CURRENT 
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Organization 
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City,  State,  Zip 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the 
New  or  Correct  Address  in  Block  6  above  and  the  Old  or  Incorrect  address  below. 
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City,  State,  Zip 


(Remove  this  sheet  along  the  perforation,  fold  as  indicated,,  staple  or  tape 
closed,  and  mail . ) 
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US  Army  Ballistic  Research  Laboratory 
ATTN:  AMXBR-OD-ST 
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